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The nonlinear theory of coUective plasmon-polariton propagation along the infinite chain of metal- 
lic nanoparticles is developed within RPA quasiclassical approach to surface plasmons in large nano- 
spheres (10 — 50 nm for radius) of Au or Ag. The wave type self- modes in the chain are determined 
and analyzed with respect to nano-sphere size and chain-separation parameters. At some regions for 
parameters the undamped modes occur. They were found on the rim of stability within the linear 
theory. The nonlinear corrections stabilize, however, diverging modes of the linear approach and 
considerably enlarge the region of undamped propagation. Nonlinearity is associated with relativis- 
tic corrections to the Lorentz friction. According to the nonlinear behavior, the region of parameters 
when undamped stable modes occur, covers the instability region of the linear theory. The remark- 
able property of fixed amplitude of the undamped collective propagating wave independently of 
initial conditions (even if they are extremely small) has been identified and described. The group 
velocity of this wave propagation is calculated. 

Keywords: plasmons, metallic nano-chain, Lorentz friction, plasmon-polariton, undamped modes 



2 



I. INTRODUCTION 



Experimental and theoretical investigations of plasmon excitations in metallic nano-crystals rapidly grew up mainly 
due to possible applications in photo- voltaics and microelectronics. A significant enhancement of absorption of incident 
light in photodiode-systems with active surfaces covered with nano-size metallic particles (of Au, Ag or Cu) with planar 
density 10*-10^°/cm^ was observed [H-Q. These findings are of practical importance for enhancement of solar cell 
efficiency, especially for developing of thin film cell technology. On the other hand, hybridized states of surface 
plasmons and photons result in plasmon-polaritons 0, Q , which are of high importance for applications in photonics 
and microelectronics 0,0, in particular, for sub-diffraction transportation of converted light energy and information 
in metallically modified structures in nano-scale [1, . 

Surface plasmons in nano-particles have been widely investigated since their classical description by Mie I 111 . Many 
particular studies, including numerical modelling of multi-electron clusters, have been carried out [H, [ij]- They 
were mostly developments of Kohn-Sham attitude in form of LDA (Local Density Approximation) or TDLDA (Time 
Dependent LDA) for small metallic clusters only [12,-16,] , up to ca 200 electrons (limited for larger clusters by numerical 
calculation constraints that grow rapidly with the number of electrons). The random phase approximation (RPA) 
was formulated (iTj for description of volume plasmons in bulk metals and utilized also for confined geometry mainly 
in a numerical or semi-numerical manner Usually, in these analyses the jellium model was assumed for 

description of positive ion background in the metal and the dynamics was addressed to the electron system only 
Such a model is preferable for clusters of simple metals, including noble metals (also transition and alkali 

metals). 

In the present paper we apply the RPA description using a semiclassical approach for a large metallic nano-sphere 
(with radius of several tens nm, and with lO'^-lO^ electrons), in an all-analytical calculus version [18]. The electron 
liquid oscillations of compressional and translational type result in excitations inside the sphere and on its surface, 
respectively. They are referred to as volume and surface plasmons. Damping of plasmons due to electron scattering 
and due to radiation losses (accounted for via the Lorentz friction force) is included. The shift of the resonance 
frequency of dipole-type surface plasmons (only such plasmons are induced by homogeneous time-dependent electric 
field), due to damping phenomena, well fits with the experimental data for various nano-sphere radii \X^- 

Collective dipole-type surface plasmon oscillations in the linear chain of metallic nano-spheres were then analyzed 
and wave- type plasmon propagation along the chain was described [20| . A coupling in the near field regime between 
oscillating dipoles in neighboring nano-spheres, together with retardation effects for energy irradiation, allowed for 
appearance of undamped propagation of plasmon waves (called plasmon-polaritons) along the chain in the experi- 
mentally realistic region of values of the separation of spheres in the chain and of the nano-sphere radii. This effect 
is of a particular significance for plasmon arranged non-dissipative and sub-diffraction transport of light converted 
energy and information along metallic chains for possible applications in nano-electronics. 

The undamped mode of plasmon-polaritons occurs, however, on the rim of stability of the linear approach. The zero 
damping rate separates the region with positive its value (corresponding to ordinary attenuation of plasmon-polaritons) 
and the region with negative damping rate (corresponding to unstable modes) . The latter exhibits unphysical behavior 
being the artefact of the linear approximation. In order to regularize the description, the nonlinear corrections must 
be thus included. The nonlinear corrections may be associated with the Lorentz friction forces. Small relativistic 
contribution to this friction has nonlinear character and quenches instable divergent modes. In the result, the insta- 
bility region of linear approach is entirely covered by the region of undamped wave propagation with the amplitude 
accommodated, however, to nonlinearity scale and independent of the initial condition (despite of its magnitude). 
This phenomenon, familiar in other nonlinear systems [21J] , seems to be of a particular significance for understanding 
of collective plasmon excitations with interesting possible applications. 



II. DAMPING OF PLASMONS IN LARGE NANO-SPHERES 

Within the RPA in semiclassical limit [31 , the solution of the dynamical equation for local density of electrons in 
a metallic nano-sphere with the radius a, can be decomposed into two parts related to the distinct domains: 

<5p(r,t) = (f^i 5"' (1) 

\Sp2{r,t), forr>a, (r ^ a-|-), ^ ^ 

corresponding to the volume and surface excitations, respectively. These two parts of local electron density fluctuations 
satisfy the equations [isj : 



dt"^ 3 m 



V^<5pi(r,t)-w^Jpi(r,t), (2) 
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and 

- m^^h{r,t) + -^^-Vj<fir,^^ (3) 
X (^pi(ri,i)e(a - ri) + 6~p2{vi,t)Q{r^ ~ a))] 5, 

where, Wp = ^^^^"^ is the bulk plasmon frequency, 8 is the Heaviside step function, 5 = 6{r — a). The analysis 
and solutions of the above equations are performed in details as presented in Ref. [isl . resulting in determination of 
plasmon self-mode spectrum, both for volume and surface modes. 

Nevertheless, this treatment did not account for plasmon attenuation. One can, however, include damping of 
plasmons in a phenomenological manner, adding an attenuation term to plasmon dynamic equations, i.e., by adding 



the term,— to the r.h.s. of both Eqs ([2|) and ([3]), taking advantage of their oscillatory form 18|. Except of 



Ait en 



To dt 

homogeneous equations ([2|) and ([3]) determining self-frequencies of plasmon modes, the dual inhomogeneous would be 
written, with explicit expression of forcing factor. This factor would be the time dependents electric field, including 
electrical component of e-m wave. For e-wave frequency in resonance with plasmons in the metalillic nanosphere, the 
wave-length (being of order of 500 nm) highly exceeds the nanosphere size (with radius 10 — 50 nm), thus the dipole 
regime is in force. For the homogeneous forcing field E(t) (which corresponds to dipole approximation satisfied for 
a ~ 10 — 50 nm, when A ~ 500 nm), only dipole surface mode can be excited and the electron response resolves to a 
single dipole type mode, described by the function Qim(t). The function Qim(t) satisfies the equation: 

a"Qi^(t) I 2 aQi^(t) I , ,2r) (f\ 

0*2 di ^^lWlm{t} 

^ [ES)5mO + V2{E^{t)5ral + Ey{t)5.m-l)\ , 

where wi = Wqi = (it is a dipole- type surface plasmon Mie frequency [HI). Only this function contributes to 
the plasmon response to the homogeneous electric field. Thus for the homogeneous forcing field, electron density 
fluctuations p^ : 

{0, r < a, 
E Qim{m^mr>a, r^a+. 
m=-l 

For plasmon oscillations given byEq. © one can calculate the corresponding dipole, 

T){t) = e J (frY5p{Y,t) = yeq(^)a^ (6) 

where, Qn(0 - /^Qxit), Qi-i(i) = ^J^Qyit), 

'^^ lx{t) and q(t) satisfies the equation (cf. Eq. (|4])), 

' q(t) = (7) 
m 
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There are various mechanisms of plasmon damping, which could be effectively accounted for via phenomenological 
oscillator type damping term. All types of scattering phenomena, including electron-electron and electron-phonon 
interactions, as well contribution of the boundary scattering effect [l^, cause significant attenuation of plasmons, 
in particular, in small metal clusters. All these contributions to damping time ratio scale as ^ and are of lowering 
significance with the radius growth. In the following subsection we argue that damping of plasmons caused by radiation 
losses scales conversely, as a'^, and for large nano-spheres this channel dominates plasmon attenuation. 



A. Lorentz friction for plasmons 



Plasmon oscillations are themselves a source of the e-m radiation. This radiation takes away the energy of plasmons 
resulting in their damping, which can be described as the Lorentz friction force reducing charge oscillations [23| . This 
damping was not included in tq in Eq. ([7]). This tq accounted only for scattering of electrons on other electrons, on 
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defects, on phonons and on nanoparticle boundary — all they lead to damping rate expressed by the simplified formula 



1 

To 



VF 

2\b 



CVf 



(8) 



where, C is the constant of unity order, a is the nano-sphere radius, is the Fermi velocity in metal. A;, is the 
electron free path in bulk (including scattering of electrons on other electrons, on impurities and on phonons p^): 
for Ag, vp = 1.4 X 10^ m/s and Af, ~ 57 nm (at room temperature); the latter term in the formula (jS]) accounts for 
scattering of electrons on the boundary of the nanoparticle, while the former one corresponds to scattering processes 
similar as in bulk. The other effects, as the so-called Landau damping (especially important in small clusters p^lzSjV 
corresponding to deca y o f plasmon for high energy particle-hole pair, are of lowering significance for nano-sphere radii 
larger than 2 — 3 nm [2J] and are completely negligible for radii larger than 10 nm. Note that the similarly lowering 
role with the radi us g rowth plays also electron liquid spill-out effect |l2l . , though it was of primary importance for 
small clusters [l^, [3 ■ 

The electron friction caused by e-m wave emission can be described as the additional electric field [23j . 



E 



L — 



(9) 



where v 



is the light velocity in the dielectric medium, and D(t) is the dipole of the nano-sphere. According to 



Eq. © we arrive at the following: 



Substituting this into Eq. ([7]), we get. 



2e Att ad^qit) 



3e?;2 3 



, 2 d 

9*2 TO dt 



q(i) 



If one rewrites the above equation (for E=0) in the form. 



9*3 ■ 



q(i) 



d_ 
dt 



To 



-q(i) 
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1 V t) / 



9*2 



thus, one notes that the zeroth order approximation (neglecting attenuation) corresponds to the equation: 

q(t) = 0. 



9*2 



(10) 



(11) 



(12) 



(13) 



In order to solve Eq. (|12p in the next step of perturbation iteration, one can substitute, in the r.h.s. of this equation, 
"!aW by -ujfq{t) (acc. to Eq. ^ 



9*2 



Therefore, if one assumes the above estimation, 
renormalized damping term: 



9jq(t) 
9*3 



,2 9q(*) 



one can include the Lorentz friction in a 



9^ 
9*2 



2 9_ 
Tdt 



q(i) = —m, 

in 



where. 



and we used for — 

To 



1 

tq 



V 



2Ab 



CvF 

2a 



(14) 



(15) 
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[221 . The renormalized damping causes a change in the shift of self- frequencies of free 
— which can be compared with the experimental observations for various nanosphere 



surface plasmons, uj[ - 
radii 

Note also, that one can verify [l^ the above calculated Lorentz friction contribution to plasmon damping by the 
estimation of the energy transfer in the far-field zone (which can be expressed by the Poynting vector) and via 
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comparison with the energy loss of plasmon osciUations. We have arrived [181 119[ at the same formula for damping 
time rate as given by Eq. ([Tl)) . The radius dependent shift of the resonance resulting due to strong irradiation- 
induced plasmon damping was verified experimentally [l9l | by measurement of light extinction in colloidal solutions of 
nanoparticles with different size (it has been done '19'] for Au, 10 — 80 nm, and Ag, 10 — 60 nm). These measurements 
clearly support the plasmon damping scaling, as described above for the far-field zone radiation losses in a dielectric 
surroundings. 

If, however, in the vicinity of the nano-sphere the another charged system is located, the situation would change. 
For instance, in the case when the nano-sphere is deposited on the semiconductor surface, the near-field coupling of 
plasmons with semiconductor band electrons must be included. 



III. ENHANCEMENT OF ENERGY TRANSFER FROM PLASMONS TO ELECTRIC RECEIVER 

LOCATED IN THE NEAR-FIELD ZONE 



Even if the derivation of plasmon dynamics equation in the form of effective harmonic oscillator equation is rigorous 
upon quantum approach of quasiclassical RPA method [18], the inclusion of plasmon attenuation of scattering type 
and of radiation losses type needs some phenomenological assumptions. They resolve themselves to extension of 
quantum RPA harmonic oscillator formulation to the damped oscillator equation form with attenuation described 
by heuristically assumed damping rates. It has been proved [l^, [2^ that radiation losses, in the case of the free 
far- field zone radiation (i.e., in the case of vacuum or dielectric surroundings of metallic nano-sphere with oscillating 
plasmons), can be accounted for as the Lorentz friction force [23l |. in the manner as described in the previous section. 
When in the near- field zone (closer than the wave length corresponding to plasmon frequency) the energy receiver (i.e., 
other system of chargesv, like semiconductor with its band system or another metallic nano-sphere as in the chain) 
is located the irradiation losses are dominated by energy transfer via this near-field zone coupling channel. Presence 
of the charged system of the energy receiver in the vicinity of e-m emitting nanosphere with plasmons, modifies a 
retarded e-m potential of the emitting system and this modifies the Lorentz friction formula, which had been derived, 
in the standard form, for the dielectric surroundings [23| . In particular, an enhancement of plasmon radiation losses in 
the case when the nanoparticles with dipole Mie surface plasmons (excited by incident external light) are deposited on 
the semiconductor surface, lies behind the observed PV efficiency growth in new generation of solar cells, metallically 
modified Jl,-^. In this case, the related attenuation rate can be also estimated by application of the Fermi golden 
rule to the semiconductor inter-band transitions induced by dipole near-field coupling with plasmons [Tsl [2^ . As it 
was proved [l8| . the resulting attenuation rate scales with nano-sphere radius, a, in different manner in comparison to 
far-field radiation, and with some correction and renormalization expressed in terms of the band system parameters 
[l^. One can expect the similar behavior in the case of the near- field coupling between nano-spheres in the chain, 
but for the sake of effectiveness of modeling one can assume that related attenuation rate has the form as that for the 
standard Lorentz friction renormalized only by some coefficient phenomenologically assumed in order to account for 
the modification of e-m potential by the receiver system presence. 



IV. NONLINEAR CORRECTIONS TO LORENTZ FRICTION FORCE 



Let us consider a metallic nano-sphere located (the center) in Rq . The electric dipole of electrons (fluctuation of 
electron density beyond the uniform distribution compensated by positive jellium) equals to. 



D(Ro,t) = e / Spir,t)rd 
Jv 



(16) 



This dipole corresponds to surface plasmons of dipole type which oscillates with Mie frequency uji — u>p/y/3e pL8| . 
where ujp is bulk plasmon frequency, e is the dielectric constant of the surrounding medium. These plasmons are not 
everlasting excitations and are damped due to scattering phenomena with the damping rate, — = tt^ + ^y^- For 
large nano-spheres the much more effective mechanism of plasmon damping are, however, irradiation energy losses, 
which for the case of irradiation to far-field zone can be expressed by the Lorentz friction [l^, [2^ . Assuming that 
electrons in the nano-sphere have positions and assuming static jellium, the dipole of the nano-sphere, D(Ro,i) = 
e^iJi-y Yi — eNeVeit), where — J2i=i ^i/^e is the mass center of the electron system. In the case of dynamics, the 
velocity of the mass center equals to, Ve = X^il^i Vj/TVe. 

On the charge eNe, located in the mass center re{t) acts a Lorentz friction force [23j . 



eA^e)' 



ds 



(17) 
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where, ds = cdt^Jl — Wg/c^, 



Up to terms of order / with respect to the niain term, one can write the electric field equivalent to the Lorentz 
friction force. 



fL _ 2 

eJVe 3 



+3^ (Ve • % 



3 d Ve ,,2 



(18) 



Next, using dimensionless variables, t' = tuji, R(t') = "iH'' -^(^') 

rj^v,(t) 



^ - R (f) = 

aJr^rft^' (*iots indicate derivatives with respect to t'), one can write out the dynamical equation in a convenient form. 
Taking into account that the dipole corresponding to surface plasmons, 

D = eNeuR., (19) 

satisfies equation of oscillatory- type, one can write it in the form (incorporating also the Lorentz friction force) , 



R (R R) + 3 R (R R)+ R (R R) 



(20) 



the terms on r.h.s. of the above equation describe the Lorentz friction including relativistic nonlinear corrections (in 
bracket) beyond the ordinary main linear term '-^R, as given by ([Q]). 

For the case when ~ ^-i^i^ ^ 1 (well fulfilled for nano-spheres with radii 10 — 50 nm, Au or Ag), one 

can apply perturbation method of solution, and in zero order perturbation assume R + R = 0. In the next step of 
perturbation one can thus substitute R= — R and R= — R in the r.h.s. of the Eq. (pOj) . 

Let us consider first a single metallic nano-sphere with dipole type surface oscillations with the dipole D. In the 
framework of the perturbation method of solution of dynamical equation of oscillatory type for the dipole, Eq. pop . 
in the first order of perturbation, attains the following form (including the damping of plasmons due to scattering 
with the rate — and due to radiation losses accounting for the linear term of Lorentz friction, while the r.h.s. of the 
equation (j2ip expresses nonlinear corrections to Lorentz friction). 



R + R 



2 

TQUll 



V3ec 



R 



/3ec 



|-f R (r- R^ + 3R (r- R^ I 



(21) 



The above nonlinear e quat ion is complicated in mathematical sense and advanced methods of solutions must be 
applied, as described in 2l|. According to the special asymptotic methods for solution of nonlinear differential 
equation (pij) . one can find the solution in the following form (R — R^), 



R{t) 



Age 



I+I7A 



o^(l-e-^) 



: COS {uJlt + 60) 



(22) 



with 



In the formulae (j22p both coefficients ^^r^- and 7 can be renormalized eventually by the mentioned above phenomeno- 
logical factor accounting for the change of^ e-m potential caused by the presence of energy receivers in near-field zone 
of plasmons on the considered metallic nano-sphere (still holding here 1 for simplicity). From the form of equation 
([22|) it follows that is always positive. Note that the scattering term, ^ = ^ + is negligible (for nano-sphere 
radius beyond 10 nm) in comparison with the linear contribution of the Lorentz friction, as it is demonstrated in Fig. 

m 
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(what is satisfied for a larger than 10 nm), 7 



1 / WpO 

™13 ivfe 



7 




Radius [nm] 

FIG. 1. Contribution to the damping rate of surface plasmon oscillations in the nano-sphere versus the nano-sphere radius, 
including the scattering attenuation (green line) and the linear Lorentz friction damping (blue line); for radii greater than ca 
10 nm the second channel dominates in overall damping (red line) 



The scale of the nonlinear corrections is given by the coefficient 7 ss 10 * y io j ■ As this coefficient is small, one 

can neglect the related contribution in the denominator for the dipole solution (|2^ . which results in ordinary linear 
solution of damped oscillations. It means that the nonlinear corrections to the Lorentz friction have no significance 
in the case of plasmon oscillations of a single nano-sphere. This situation changes, however, considerably in the case 
of collective plasmon excitation propagating along the metallic nano-chain, as it will be described in the following 
paragraph. 



V. COLLECTIVE PLASMON WAVE-TYPE PROPAGATION ALONG THE NANO-CHAIN IN THE 

NONLINEAR REGIME 



In the case of the metallic nano-chain one has to take into account the mutual affecting of nano-spheres in the 
chain. Assuming that in the sphere located in the point r we deal with the dipole D, then in the other place rg (ro is 
fixed to the end of r) the dipole type electric field attains the form as follows (including electro-magnetic retardation). 



E(r,ro,i) = ^{3no (no- D (r, t - ni)) 



(23) 



This allows for writing out the dynamical equation for plasmon oscillations at each nano-sphere of the chain, which 
can be numbered by integer I (d will denote the separation between nano-spheres in the chain; vectors r and rg are 
coUinear, if the origin is associated with one of nano-spheres in the chain). The first term of the right-hand-side 
in the following formula (|24l) describes the dipole type coupling between nano-spheres [13] and the other two terms 
correspond to contribution due to plasmon attenuation (in the latter term the Lorentz friction caused electric field 
accounts also for nonlinear corrections). The index a enumerates polarizations, longitudinal and transversal ones with 
respect to the chain orientation. 



Ra + Ra{ld,t) — ^ 

m— — C30,m^/ 

--l-R^[ld,t) + ^E^{ld,t), 



(24) 



I — Id X v , , 

where, '^a ~ { 2 a — z introduced to distinguish both polarizations. The summation in the first term of the 

r.h.s. of the equation (j24p can be explicitly performed in the manner as presented in [20j . because it is the same as for 
the linear theory formulation. Similarly as in the linear theory framework, one can change to the quasi-momentum 
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FIG. 2. Regions for negative value of damping rates for plasmon-polaritons in the chain (1 for longitudinal polarization modes 
and 2 for transversal one, in the nonlinear formulation framework; for the linear theory, red lines gives the position of vanishing 
damping rate for longitudinal modes of plasmon-polaritons and blue lines the same for transversal modes 



picture, taking advantage of the chain periodicity (in analogy to Bloch states in crystals with the reciprocal lattice of 
quasi- momentum), i.e.. 



< k < ^,Ra {k) sin {toji + /3) 



2. 6 , «^ (25) 



Thus the equation (I24p can be rewritten in the following form (the Lorentz friction term was represented similarly as 
in equation (|211) '). 



where, ^ 1 - 2crQ^ cos (kd) cos (^) and 

- U^) ^ ( (^f ~ (kd- nf + 4) , 

(27) 



In the above formulae the remarkable property is linked with the expressions for the attenuation rate for both 
polarizations. Two last expressions below equation (j26p give these damping rates explicitly and one can notice that 
they could change their signs depending on values for d, a and k. In Fig. ([2]) the regions of negative value for damping 
rates are marked (for both polarizations). 

Applying the same methods for solution of the nonlinear equation ()26p as in the former paragraph, using the 
asymptotic methods (2]| , one can find the corresponding solutions for both regions with positive and negative damping 
rate, respectively. 

For the positive damping rate, ;^r^ > 0, 

^ 

i?a {k, t) = — ^sS^^S cos [uJo^t + Go) , 

Wl+7.A^„o(l--"^) (28) 
Ra {k,t) -J'(t_>oo) 0, 
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where, 7^ = ItqWiI (^^)^ j (i'^a ~ l)- We note from the form of Eq. (|28p that this is a damped mode, vanishing at 
longer time scale. 

Nevertheless, for negative damping rate, -:f^ < 0, the solution has a different form. 



Ra (fc, t) = — COS {uJct + 60) 



(29) 



Ra (fc, t) ^(t_,oo) -7^^ COS (tWa + Oq) . 



This solution is stable. It corresponds to an undamped mode which stabilizes on the fixed amplitude at longer time 
scale, independent of initial conditions expressed by Aa ■ 

The corresponding dipole oscillations attain in the latter case the form of 'planar' waves propagating along the 
chain. 



_ eN,a 1 

J-^n — 



/7a 



From the above formulae it follows that for positive attenuation rate we deal with ordinary damped plasmon- 
polariton propagation, not strongly modified in comparison to linear theory (due to small value of the factor 7). 
Nevertheless, in the case of negative damping rate the solution behaves differently-on longer time scale this solution 
stabilizes on the constant amplitude independently of initial conditions. This remarkable property characterizes 
undamped propagation of plasmon-polariton along the chain. If one turns back to dipole explicit form (|19p. then 
typical 'planar' wave formula with constant amplitude describes this undamped mode, as written in the equation (|30l) . 
The region of negative damping correspond thus, within the nonlinear approach, to undamped modes with the fixed 
amplitude. Let us note that the same region was linked with instability of the linear theory (which was, however, the 
artefact of the linear approach) . 

Finally, one can calculate the group velocity of the undamped plasmon-polariton mode, in the following form, 

diUa , a^a^ sin {kd) COS {^) 
Va^ — ^UJld . ^ ^=. (31) 

dyi-2(7„f^cos(W)cos(^^) 

From this formula it follows that the group velocity of the unadamped wave type collective plasmon excitation (called 
plasmon-polariton) may attain different values depending on a, d and k. 

Indicated above undamped mode of propagation of collective surface plasmons seems to match with experimentally 
observed log range propagation of plasmon excitations along the metallic nano-chain 8, 10, 22, 26, 27]. The constant 
and fixed value of the amplitude for these oscillations ([25]) are independent of initial conditions, which means that 
these excitations will be present in the system even if are excited by arbitrary small fluctuations. Thus one can 
conclude that they are self-exciting modes which always present in the system provided that radii of spheres and their 
separation in the chain have values for which at least one of the attenuation rates (j27p is negative. 



VI. CONCLUSIONS 



We have demonstrated the practical utilization of RPA semiclassical description of plasmon oscillations in metallic 
nano-spheres. The oscillatory form of dynamics both for volume and surface plasmons, rigorously described upon 
the RPA semiclassical limit fits well with the large nano-sphere case, when the nano-sphere radius is greater than 
10 nm and lower than 60 nm, (for Au, Ag or Cu material), what is confirmed by experimental observations, on the 
other hand. The most important property of plasmons on such large nano-spheres is the very strong e-m irradiation 
caused by these excitations, which results in quick damping of oscillations. The attenuation effects for plasmons were 
not, however, included into the quantum RPA model. Nevertheless, they could be included by a phanomenological 
manner, taking advantage of the oscillatory form of dynamical equations. Some information on plasmon damping 
can be taken from microscopic analyzes of small metallic clusters (especially made by LDA and TDLDA methods of 
numerical simulations employing Kohn-Sham equation). For larger nano-spheres, these effects, mainly of scattering 
type (also Landau damping), are, however, not specially important as diminishing with radius growth, as i. 

The irradiation effects overwhelming the energy losses in the case of large nano-spheres can be grasped in terms 
of the Lorentz friction, which reduces the charge movement. This approach has been analyzed in the present paper. 
Two distinct situations were indicated, the first one-of the free radiation to far-field zone in dielectric (or vacuum) 
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surroundings of single nanoparticle and the second one, when in the near-field zone of plasmons on the nano-sphere, 
an additional charged system is located. 

This additional charge system acting as the e-m energy receiver, strongly modifies the e-m potentials of the source 
and in this way modifies energy emission in comparison to the free emission in vacuum or in dielectric surroundings. In 
particular, the Lorentz friction is modified in the case of energy receiver presence in the near field zone of plasmonms, 
in comparison to simple free emission to the far-field zone. The e-m energy receiver located close to emitting nano- 
sphere, could be semiconductor (as in the case of metallically modified solar cells) or other metallic nano-spheres (as in 
the case of metallic nano-chain). The latter situation has been analyzed in this paper. We have shown previously [20j 
that along the infinite nano-chain the collective plasmon-polaritons can propagate (being collective surface plasmons 
coupled by e-m field in near-field zone), which at certain values of nano-sphere radii and separation in the chain, 
appear as undamped modes. Simultaneously, the instability regions of linear theory of plasmon-polariton dynamics 
occur, which shows that the nonlinear corrections must be included. 

In this paper we have developed the nonlinear theory of collective plasmon-polariton dynamics along the chain, 
including nonlinear corrections to Lorentz friction force. Even though the related nonlinearity is small, it suffices 
to regularize the instable linear approach. As the most important observation, we noted the presence of undamped 
excitations (instead of those instable within the linear approach), which have fixed amplitude independently how 
small or large the initial conditions were. This excitations, typical for nonlinear systems, would have some practical 
significance, e.g., to enhance sensitivity of antennas with coverings by plasmon nano-systems offering self- induced 
collective plasmon-polaritons in wide range of frequencies, which would be excited by even very small signal (the 
energy to attain the stable level of plasmon-polariton amplitude would be supplied, in this case, by an external 
auxiliary supply). 



[1] K. Okamoto, I. Niki, A. Scherer, Y. Narukawa, Y. Kawakami, Appl. Phys. Lett. 87, 071102 (2005) 
[2] S. Pillai, K.R. Catchpole, T. Trupke, G. Zhang, J. Zhao, G. M.A, Appl. Phys. Lett. 88, 161102 (2006) 
[3] D.M. Schaadt, B. Feng, E.T. Yu, Appl. Phys. Lett. 86, 063106 (2005) 
[4] S.P. Sundararajan, N.K. Grandy, N. Mirin, N.J. Halas, Nano Lett. 8, 624 (2008) 

[5] M. Westphalen, U. Kreibig, J. Rostalski, H. Liith, D. Meissner, Sol. Energy Mater. Sol. Cells 61, 97 (2003) 

[6] A.J. Morfa, K.L. Rowlen, T.H. Reilly, M.J. Romero, J. Lagemaat, Appl. Phys. Lett. 92, 013504 (2008) 

[7] S.A. Maier, Plasmonics: Fundamentals and Applications (Springer, Berlin, 2007) 

[8] A.V. Zayats, I.I. Smolyaninov, A.A. Maradudin, Phys. Rep. 408, 131 (2005) 

[9] F.J.G. de Abajo, Rev. Mod. Phys. 82, 209 (2010) 
[10] S.A. Maier, P.G. Kik, H.A. Atwater, Phys. Rev. B 67, 205402 (2003) 
[11] G. Mie, Ann. Phys. 25, 337 (1908) 
[12] M. Brack, Rev. of Mod. Phys. 65, 667 (1993) 
[13] M. Brack, Phys. Rev. B 39, 3533 (1989) 
[14] W. Ekardt, Phys. Rev. B 31, 6360 (1985) 
[15] W. Ekardt, Phys. Rev. B 33, 8803 (1986) 
[16] V.V. Kresin, Phys. Rep. 220, 1 (1992) 

[17] D. Pines, Elementary Excitations in Solids (ABP Perseus Books, Massachusetts, 1999) 

[18] J. Jacak, J. Krasnyj, W. Jacak, R. Gonczarek, A.Chepok, L. Jacak, Phys. Rev. B 82, 035418 (2010) 

[19] W. Jacak, J. Krasnyj, J. Jacak, R. Gonczarek, A. Chepok, L. Jacak, D. Hu, D. Schaadt, J. Appl. Phys. 107, 124317 (2010) 
[20] W. Jacak, J. Krasnyj, J. Jacak, A. Chepok, L. Jacak, W. Donderowicz, D. Hu, D. Schaadt, J. Appl. Phys. 108, 084304 
(2010) 

[21] N.N. Bogolubov, J. A. Mitropolskyj, Asymptotical methods in theory of nonlinear oscillations (Nauka, Moscow, 2005) 
[22] M.L. Brongersma, J.W. Hartman, H.A. Atwater, Phys. Rev. B 62, R16356 (2000) 
[23] L.D. Landau, E.M. Lifshitz, Field Theory (Nauka, Moscow, 1973) 

[24] C. Yannouleas, R.A. Broglia, M. Brack, P.P. Bortignon, Phys. Rev. Lett. 63, 255 (1989) 

[25] W. Jacak, J. Krasnyj, J. Jacak, W. Donderowicz, L. Jacak, J. Phys. D: Appl. Phys. 44, 055301 (2011) 

[26] V.A. Markel, A.K. Sarychev, Phys. Rev. B 75, 085426 (2007) 

[27] W.L. Barnes, A. Dereux, T.W. Ebbesen, Nature 424, 824 (2003) 



